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Abstract 

An exactly solvable model suitable for the description of single 
and double-beta decay processes of the Fermi-type is introduced. The 
model is equivalent to the exact shell-model treatment of protons and 
neutrons in a single j-shell. Exact eigenvalues and eigenvectors are 
compared to those corresponding to the hamiltonian in the quasipar- 
ticle basis (qp) and with the results of both the standard quasiparti- 
cle random phase approximation (QRPA) and the renormalized one 
(RQRPA). The role of the scattering term of the quasiparticle hamilto- 
nian is analyzed. The presence of an exact eigenstate with zero energy 
is shown to be related to the collapse of the QRPA. The RQRPA and 
the qp solutions do not include this zero-energy eigenvalue in their 
spectra, probably due to spurious correlations. The meaning of this 
result in terms of symmetries is presented. 
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1 Introduction 

In the last years the study of the Quasiparticle Random Phase Approxima- 
tion (QRPA) and its extensions, like the Renormalized Quasiparticle Random 
Phase Approximation (RQRPA), has received renewed attention. The goal 
was to improve substantially the reliability of the QRPA description of nu- 
clear double beta decay transitions and, at the same time, to enhance the 
predictive power of the theory in an unambiguous way 

The predictive power of the QRPA, mostly in dealing with the calculation 
of the matrix elements for ground state to ground state two-neutrino double- 
beta decay transitions ( /3/3 2i/ ), is questionable since these amplitudes are 
extremely sensitive to details of the nuclear two-body interaction |1|, |2], [3], |]. 

The inclusion of renormalized particle-particle correlations in the QRPA 
matrix amounts to a drastic suppression of the fifav -matrix elements. How- 
ever, for some critical values of the model parameters; i.e: the renormalized 
two-body interactions, the otherwise purely real QRPA eigenvalue problem 
becomes complex. As a consequence of it the standard properties of the 
QRPA metric and conservation rules are severely downplayed by the ap- 
pearance of strong ground-state correlations which jeopardize the stability 
of the theory. The most notorious example of this behaviour, of the QRPA 
approach, is the calculation of the fifiiv decay of 100 Mo |2|, ||, ||, [/J. 

The renormalized version of the QRPA (RQRPA) || [|, which includes 
some corrections beyond the quasiboson approximation, has been recently 
reformulated [JTOj and applied to the Pfov decay problem Contrary to 



the QRPA, the RQRPA does not collapse for any value of the residual two- 
body interaction. Based on its properties, the RQRPA was presented as a 
cure for the instabilities of the QRPA and it was applied to calculations of 
the (3(32v decay of 100 Mo f[T |. Similar studies have been performed in the 



framework of the RQRPA and with the inclusion of proton-neutron pairing 
correlations in symmetry breaking Hamiltonians. ||12|| . 

In a recent paper |H| we have shown that the RQRPA violates the Ikeda 
sum rule and that this violation is indeed present in many extensions of the 
QRPA. The study was based on the schematic proton-neutron Lipkin model. 

In a subsequent letter III we have introduced an exactly solvable model 
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for the description of single- and double-beta decay Fermi-type transitions. 
This model is equivalent to a single-j shell model for protons and neutrons. 
The appearance of an eigenvalue at zero energy, in the exact spectrum, was 
found. Moreover, it has been shown that the presence of this zero-energy 
eigenvalue should be associated to the collapse of the QRPA. It was shown 
that the RQRPA does not include this zero-energy mode in its spectrum. It 
was also shown that the absence of this zero-energy state, in the RQRPA, 
leads to finite but spurious results for the transition amplitudes near the 
point of collapse of the QRPA. 

In the present paper we discuss the details of the exactly solvable model 



of ||14j| . The algebraic techniques needed to evaluate matrix elements of the 
relevant operators, in the SO (5) group-representation, are described in detail. 
Exact eigenvalues and eigenvectors are compared with those corresponding to 
the quasiparticle version of the hamiltonian (qp) and with the ones obtained 
with the QRPA and RQRPA. The role of the correlations induced by the 
scattering term H31 of the qp-hamiltonian and the effects on the number 
of quasiparticles in the ground state are analyzed. The presence of a zero 
excitation energy state in the spectrum corresponding to the exact solution 
of the model hamiltonian is discussed. As said before it will be shown that 
the RQRPA and the qp solutions, do not display the same feature, most 
likely due to the presence of spurious states caused by the mixing of orders, 
of the relevant interaction terms, in the expansion procedure. 

The structure of the paper is the following: the model and its solutions are 
presented in Section 2, the quasiparticle version of the hamiltonian, its linear 
representation in terms of pairs of unlike (proton-neutron) quasiparticle-pairs 
and its properties are introduced in Section 3. The QRPA and RQRPA treat- 
ments of the hamiltonian are discussed in Section 4. The matrix elements 
of double-beta-decay transitions, calculated in the framework of the different 
approximations introduced in the previous sections, are given in Section 5. 
Conclusions are drawn in Section 6. The SO(5) algebra, representations and 
reduced matrix elements used in the calculations are given in detail in the 
Appendices A, B and C, respectively. 
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2 The model 



The model hamiltonian, which includes a single particle term, a pairing term 
for protons and neutrons and a schematic charge- dependent residual inter- 
action with both particle-hole and particle-particle channels, has been intro- 
duced in refs. flU^, [16|, |17j and it is given by 



H = e v M v - G V S% + e n M n - G n S\S n + 2 X /T • (3 + - 2kP~ • P+ (1) 
with 



M = E o) mi a mt , Sj = E «U a k/ 2 

rrij mi 

= E a m„ a "in ) -P = E a m„ G m 



p, n 



(2) 



m v =m n r m„=— m„ 



aj, = a] pmp being the particle creation operator and ajj = (— ^Y p ~ mp a ] v ~m p its 
time reversal. The parameters x and k play the role of the renormalization 
factors g p h and g pp introduced in the literature |], |[ |J. 

It has been shown in a series of papers |T^, |1^, |T^] that this hamiltonian, 
when treated in the framework of the QRPA, reproduces fairly well the re- 
sults obtained with a realistic G-Matrix constructed from the Bonn-OBEP 
potential, both for single- and double-beta decay transitions. These results 
can be taken as an indication about the correlations induced by the interac- 
tions in (1), which are obviously specific to the relevant degrees of freedom 
of the problem. In other words, if the relatively simple schematic force (1) 
can approximately described the correlations induced by a more realistic in- 
teraction it certainly means that it is able to pick-up the bulk of the physics 
involved in the transitions. 

In a single-one-shell limit, for the model space (j p = j n = j) and for 
monopole (J = 0) excitations the hamiltonian (1) can be solved exactly. In 
spite of the fact that the solutions obtained in this restricted model space 
cannot be related to actual nuclear states, the excitation energies, single- 
and double-beta decay transition amplitudes and ground state correlations 
depend on the particle-particle strength parameter k in the same way as 
they do in realistic calculations with many single particle levels and with 
more realistic interactions, as we shall show later on. Physically, the beta 
decay transitions between J 71 = + states correspond to transitions of the 
Fermi type. However, the study of the model and the identification of its 
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relevant degrees of freedom, instead of the comparison of observables, is the 
main aspect of the present work. We shall obtain the eigenstates of (1), by 
using different approximations, in order to built-up a comprehensive view 
about the validity of them and their predictive power. 

The hamiltonian ([]]) can be expressed in terms of the generators of an 
SO (5) algebra |20], |21], p2] . The Hilbert space is constructed by using the 
eigenstates of the particle-number operator Af = Af p + Af n , the isospin T 
and its projection T z = (Af n — Af p )/2. The raising and lowering isospin 
operators are defined as (3 = T , where T~\n) = \p). With them we can 
construct the isospin scalar T 2 = \{T~T + + T +T~) + T z and the second 
order SO (5) Casimir (see Appendix A) 

SlS n + S}£ p + Vp = ^(3 - ^ + 2fi) - ^(T + 1) (3) 

with tt = (2j + l)/2. 

The Hamiltonian (|l]) can be expressed in terms of the above mentioned 
operators. Hereafter we will use G p = G n = G for simplicity. In terms of 
these generators the Hamiltonian (1) reads 

H= [e p + e n -i(3 + 2n-^)(G + 2K)]f + 
[e„ - e p - 2x(T z - 1)]% + 

[2 X + f + |«]T(T + 1) + ( 4 ) 



3 1 3 

tt(4:K - G) 



r a t at iJ=o,r=ir aal J=o,r=i| 



r=2 
r z =o 

In writing the creation and annihilation operators (a^ , a) we have omitted 
unnecessary subindexes since the coupling to total angular momentum J and 
isospin T, represented as [aW] J,r , is understood. 

Hamiltonian (|4j) is diagonal in the A/", T, T z basis if G = 4/t. It can be 
reduced to an isospin scalar if its parameters are selected as 

e P = e n , x = 0, G = 4k. (5) 

If 4k ^ G the hamiltonian (JXj) is not diagonal in this basis. The hamil- 
tonian mixes states with different isospin T while its eigenstates still have 
definite N and T z . The dynamical breaking of the isospin symmetry is an 
essential aspect of the model which is directly related to the nuclear struc- 
ture mechanism responsible for the suppression of the matrix elements for 
double-beta-decay transitions. 
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2.1 The diagonal case G = 4k, 

The solution of (1) in the basis | Af,T,T z >, in the case G = 4n, gives a 
state, the isobaric analog state (IAS) at the energy 

E IAS = E(N,T,T z = T-\)-E(N,T,T z = T) 
= e p -e n + 4 X (T-l) 

Considering a double Fermi-transition, the energy available for the decay 
is given by 

En =E(Af,T,T z )-E(Af,T,T z ) m 
= 2(e n -e p )+8x + G(2T-l) , if T > 2 . { ' > 

where T z = T z — 2. The above expression shows clearly the role of the 
particle-hole strength parameter x- It determines the excitation energy of 
the IAS, which depends not only upon the proton-neutron energy shift due to 
the nuclear Coulomb field but also upon x an d T. The same dependence is 
shown by the Q value (eq.(7)). In analogy with the situation found in realistic 
calculations its value can be determined by a fit to the experimental value 
of the IAS energy (or to the Gamow- Teller Giant Resonance for the case of 
J 71 = 1 + spin-isospin-dependent excitations). 

The P decay operators for single Fermi transitions, T ± , do not change 
the total isospin neither the total particle number of the state upon which 
they act. Only the isospin projection of the state is changed in steps of one 
unit, namely: 

T± \MTT Z ) = ^J(T±T Z + 1)(T T T Z ) \MTT Z ± 1) (8) 

2.2 The spectrum 

For the numerical examples we have selected N n > N p and a large value 
of j to simulate the realistic situation found in medium- and heavy-mass 
nuclei. To perform the calculations we have adopted the following two sets 
of parameters: 
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set I 



3 = 9/2, 

e p = 0.96MeV, 

G p = G n = OAMeV, 



X 



e 



N 



n 



10, 

: O.OMeV, 

or O.OAMeV , 



< k < 0.2 



< T z < 4, 



and 
set II 



(9) 



3 = 19/2, 

e p = O.Q9MeV, 

G p = G n = 0.2MeV, 



X 



N 



e, 



n 



20, 

: O.OMeV, 

or 0.02hMeV , 



< k < 0.1 



1 < < 5, 



The dependence of the spectrum and transition matrix elements on the 
parameters x an d k is analyzed in the following paragraphs. 



The complete set of + states, belonging to different isotopes, is shown in 
Fig. la and Fig. 2a, for G — 4k and x — 0., as a function of the number of 
protons (Z). The states are labeled by the isospin quantum numbers (T,T Z ). 
Ground states are shown by thicker lines. As shown in these figures the 
structure of the mass parabola is qualitatively reproduced. 

The upper insert, case a) of each figure, shows the full spectrum corre- 
sponding to x — 0- The lower one, case b), shows the results corresponding 
to x = O.OhMeV (Fig Lb) and x = 0.025MeV in (Fig 2.b). Obviously the 
particle-hole channel of the residual interaction stretches the spectra of all 
isotopes. As mentioned above, it increases the energy of the IAS. 

Beta decay transitions of the Fermi type, mediated by the action of the 
operator f3~ = t _ , are allowed between states belonging to the same isospin 
multiplet. The energy of each member of a given multiplet increases linearly 
with Z. 

In this example the + states belonging to each odd-odd-mass nuclei ( N- 
1, Z+l, A) are the IAS constructed from the + states of the even-even-mass 
nuclei with (N, Z, A) nucleons. Thus, Fermi transitions between them are 
allowed. 

Since the isospin of the ground state of each of the even-even-mass nuclei 
differs, for different isotopes, Fermi-double-beta-decay transitions connecting 
them are forbidden in this diagonal limit G = 4k. 



Fig. 1 
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2.3 Exact solutions 



The Hamiltonian (|I]) has a T = 2 tensorial component which mixes states 
with different isospin, while particle number and isospin projection remain as 
good quantum numbers. The diagonalization of (1) is performed in the basis 
of states described in the Appendix B. The corresponding reduced matrix 
elements are given in the Appendix C. The eigenstates are written as: 



\NT z a) = J2 C AfTT z \NTT z ) (10) 

T 

Fig. 3 

The energy of the ground-state (0+J and of the first-excited state (0+), 
as a function of the ratio Ak/G for the set of parameters j = 9/2, M n = 
6, M p = 4, x = are shown if Fig. 3a. The results of Fig.3b have been 
obtained with the set of parameters given by j = 19/2, Af n = 12, Af p = 8 
and x — 0- 

The most characteristic feature of the results is the barely avoided crossing 
of levels, due to the repulsive nature of the effective residual interaction 
between them. Although a complete level-crossing is not obtained in this 
model, in the neighbourhood of the value 4k/ G ~ la major structural change 
in the wave functions will develope. In the case of a complete crossing of 
levels the ground and the first excited state will interchange their quantum 
numbers thus given raise to a permanently deformed (in the sense of the 
isospin dominance) situation. 

This behaviour is by no means a surprise since it is similar to that found 



in pairing plus quadrupole systems [23[. In this case, if the quadrupole 



quadrupole interaction is strong enough, the system becomes permanently 
deformed, in the sense of the angular momentum and spatial rotations. The 
analogy between this and the present case ( isospin degree of freedom) can be 
drawn from the study of PIL |25| where the "pairing plus monopole" model, 



which is a two- level model exactly solvable using the SO (5) algebra, was 
used to analyzed the spherical and the deformed regime of the solutions of 
the multipole-multipole interaction. 

Figs. 4,5 

The full-thin line of Fig. 4a (4b) represents the excitation energy E exc of 
the lowest + state belonging to the double-odd-mass nucleus (jV n = 7, Af p = 
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3) with respect to the parent even-even-mass nucleus (Af n = 8,Af p = 2) as 
a function of the ratio 4k /G for j = 9/2 and x = (0.04). It is clear that 
when 4k/ G ~ 1.6(1.8) attractive proton-neutron correlation dominates over 
proton-proton and neutron-neutron pairing correlations and the excitation 
energy goes to zero. Similar results are depicted in Figs. 5a and 5b, corre- 
sponding to the excitation energy E exc of the lowest + state in the odd-odd 
mass nucleus (An = 13, Af p = 7), also measured from the ground state of 
the parent even-even nucleus with (Af n = 14, A/^ = 6) , for j = 19/2 and 
X = (0.025). In the case of Fig. 5 the excitation energy goes to zero when 
4k/Gk 1.3. 

The vanishing of the energy of the first excited state, and the subsequent 
inversion of levels (or negative excitation energies) would indicate that the 
double-odd nucleus becomes more bound than their even-even neighbours, 
contradicting the main evidence for the dominance of like-nucleons pairing 
in medium- and heavy-mass nuclei. It would also completely suppress the 
double beta decay because the single beta decay from each "side" of the 
double-odd nucleus would be allowed. 

These result simply emphasizes the fact that the Hamiltonian ([I]) will 
not be the adequate one when attractive proton-neutron interactions are 
too large. In a realistic situation, obviously, the true Hamiltonian includes 
other degrees of freedom, like quadrupole-quadrupole interactions, and per- 
manent deformations of the single-particle mean-field can also be present. 
These additional degrees of freedom will prevent the complete crossing of 
levels which, of course, is not observed. However, in many cases the exper- 
imentally observed energy-shift of double-odd-mass nuclei, respect to their 
double-even-mass neighbours is very small. This finding reinforces the notion 
of an underlying dynamical-symmetry-restoration-effect. 



3 The Hamiltonian in the quasiparticle (qp) 
basis 

By performing the transformation of the particle creation and annihilation 
operators of the Hamiltonian (1) to the quasiparticle representation J2| 



i.e. 



by using the Bogolyubov transformations for protons and neutrons, we have 
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obtained the Hamiltonian 



H = (e p - \ P )N P + (e n - A n )iV n + \ X A*A + A 2 (A I "A I " + AA) 

-\ 3 (A^B + B*A) - A 4 (A t 5 t + BA) + X 5 B^B + A 6 ( J B t 5 t + 55) 

(11) 

where e p , e n are the quasiparticle energies, A p , A n the chemical potentials and 



A* = \al (g) at 



j=o 

M=0 



B^ = \al (g) a f 



j=o 

M=0 



t a- 



i=p,n 



Ai = An [x{u 2 p v 2 n + v 2 p u 2 n ) - K{u 2 p u 2 n + v 2 p v 2 n )\ , A 2 = 4fi(x + n)u p v p u n v n 
A3 = 4f2(x + K)u n v n (u 2 p - uj) , A 4 = 4fi(x + K)u p v p (u 2 n - v 2 ) , 



Af; — — Ao 



(12) 

The operators A* (A), which create (annihilate) a pair of unlike (pro- 
ton and neutron)-quasiparticles, together with their counterparts for pairs 
of identical quasiparticles and B, B\ N p , N n are the generators of the SO(5) 
algebra Pfl . 

The quasiparticle energies 



GQ/2 . 



(13) 



and the occupation probabilities 

AC, 



2j + l 



■A4 
2j + l 



(14) 



are determined from the gap equations and particle-number conservation 
condition p6| . The occupation probabilities can also be defined in terms of 
the single-particle and quasiparticle energy, namely: 



1 _ e, - Gvf - Xj 
2 [ u j ' 



p, n 



(15) 



From this equation and from Eq. ([13]) the chemical potentials can be 
expressed as 



A; 



Gtt 



e,- - 



+ Gv 



p, n 



(16) 
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The excitation energy, E^ xc , of a state |0a) belonging to the spectrum of a 
double-odd mass nucleus, with J\f p + 1 protons and Af n — 1 neutrons, respect 
to the ground state of the even-even neighbour with Af p ,J\f n , can be easily 
calculated if blocking is considered, i.e. when v p , v n are calculated for the 
even-even and odd-odd nuclei separately. These excitation energies are given 
by 

E x exc = E(\,J\f p + 1, M n - 1) - E{g.s.,M n ,M p ) + X P -\n (17) 



In the following we shall always refer to Eq. ( Jl7|) as a suitable approxima- 
tion for the excitation energies. In the present calculation we have selected 
e p — e n is such a way that 



G Q — 1 

\-K = e p -e n - -{M n - K)—£f- = > ( 18 ) 



which implies 



G Vt — 1 

Alternatively, one can compute the occupation amplitudes v p , v n always 
for the even-even nucleus, without including blocking. The effect of blocking 
on the unperturbed excitation energies, with k = \ = 0, can be ignored if the 
single-particle energy difference between protons and neutrons is modified to 
the value 

e P - e n = ^(Af n -J\f p — l)^y^ • (20) 

The linearized version of the Hamiltonian (0) is obtained by keeping only 
the first line of Eq. ([11]). This is equivalent to neglect terms proportional to 
B and (the so-called scattering terms). The solutions of this truncated 



Hamiltonian have been discussed in a previous paper JT3[ . 

Finding the eigenvalues and eigenvectors of Hamiltonian (11) requires 
the use of the same algebraic techniques involved in solving the original 
Hamiltonian. However, the complexity of the problem increases severely, due 
to the fact that neither the quasiparticle number or the quasiparticle isospin 
projection (or equivalently the number of proton and neutron quasiparticles) 
are good quantum numbers. It implies that the dimension of the basis will 
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increase by two orders of magnitude. Additional reduced matrix elements 
would then be needed to diagonalize the Hamiltonian (|TTD . The analytic 
expressions of these matrix elements are given in Appendix C. 

There is a remaining symmetry in Hamiltonian (PD, since states with 
even number of proton- and neutron-quasiparticles are not connected with 
states having an odd number of them. Due to this fact it is possible to 
diagonalize separately these two cases. 

Particle number is not a good quantum number, obviously, because it 
is broken spontaneously by the Bogolyubov transformation. Thus, zero- 
quasiparticle states belonging the even-even-mass nucleus have good average 
number of protons and neutrons, the condition used to determine v p ,v n , while 
states with a non- vanishing number of quasiparticles show strong fluctuations 
in the particle number. Fluctuations in the particle number can induce, 
naturally, important effects on the observables. Moreover, the admixture 
of several quasiparticle-configurations in a given state, induced by residual 
particle-particle interactions, can also strongly influence the behaviour of 
the observables. An example of this effect is given in [fQ5| , concerning the 
violation of the Ikeda Sum Rule produced by large values of the particle- 
particle strength k. 

The spectrum of the qp-Hamiltonian (|TTD is shown in Figures 4 and 5. 
The curves shown by small-dotted-lines, in Figs. 4. a), 4.b), 5. a), 5.b), display 
the dependence of the excitation energy for the qp-hamiltonian (|ll]) upon the 
ratio k/G. The results of this qp-approximation closely follow the exact ones 
up to the point where they become negatives (4k/ G ~ 1.4 — 1.8 in the 
different cases). From this point on they vanish, rather than taking negative 
values, instead. The excitation energies for the linearized hamiltonian H22 + 
ifo4 are shown as thick lines in these figures. We can see that the linearized 
hamiltonian is able to reproduce qualitatively the behaviour of the full-qp 
one, but in general it overestimates the values of the excitation energies. 

As it is mentioned above, the results shown in Figures 4 and 5 have been 
obtained both with the complete qp-hamiltonian and with the truncated 
hamiltonian which includes only the product of pair-creation and annihilation- 
operators. In |13| the relevance of the scattering terms in (|Tl~D was pointed 



out. From the present results it can be seen that the inclusion of these terms 
is indeed important if one looks after a better description of the qp-excitation 
energies, up to the point where the exact excitation energies become nega- 
tive. For larger values of k even the eigenstates of the complete hamiltonian 
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fail to describe negative excitation energies. This is a clear indication that 
other effects can play an important role, like, i.e; effects associated to the 
appearance of spurious states. This can be quantitatively illustrated by the 
following. There are four exact eigenstates for j = 19/2, Af n = 13, Af p = 7, 
as can be seen in Fig. 2. a), while the spectrum of the qp-hamiltonian fllTD has 
220 eigenstates. It is well known that states with M n = 14±iV n , M p = 6±N P , 
where N p and N n are the number of quasiparticle protons and neutrons, re- 
spectively, are mixed with two-(p-n)-quasiparticle states in the odd-odd nu- 
cleus and provide a large number of states belonging to other nuclei. When 
4k/G <ti 1 the spurious states remain largely un-mixed with the lower energy 
two-qp state. But when 4k,/ G ~ 1 the mixing becomes important. This fact 
up-grades the relevance of particle-number violation effects in dealing with 
this case. 

The full qp-treatment represents the best possible extension of the quasi- 
boson approximation, without performing a particle-number projection, in a 
single-j shell. It goes beyond any second extended RPA p7j and it includes 



explicitly all number of proton- and neutron-quasiparticles (N p and N n ) in 
the eigenstates. 

To analize the effects associated to the number of quasiparticles in the 
ground state of double-even nuclei, and particularly the effects associated 
to the number of quasiprotons, we have calculated the average number of 
quasi-protons using the expression 

(o A |i\go A > = ]T \c x NTT f(N/2 + T z ) . (21) 

NTT Z 

A similar expression holds for the average neutron-quasiparticle-number. 
Fig 6, 7 

In Figs. 6. a) and 6.b) the average number of proton-quasiparticles in 
the ground state of the even-even nucleus with Af p = 6, Af n = 14 is shown 
as a function of 4k/G, for j = 19/2, x = and 0.04. Figs. 7.a) and 
7b) show the number of proton-quasiparticles for Af p = 8, Af n = 12, j = 
19/2, x — an d 0.025. The dashed-lines represent the results corresponding 
to the full qp-hamiltonian case while the large dots refer to the linearized 
H22 + H04 version of it. The difference between both approximations is 
evident. Using the linearized hamiltonian the states are composed only by 
proton-neutron-quasiparticle pairs ||13|| , while the presence of the scattering 



terms introduces also like-(p-p and n-n)-quasiparticle pairs. The presence of 
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these pairs, which for 4k/ G ~ 1 play a crucial role, increases notably the 
number of quasiparticles and yields excitation energies closer to the exact 
ones. 

The average quasiparticle number shows a saturation in the full-qp case 
for 4k/G ~ 1.8. At this value of the residual pn-interaction the ground state 
is far-away for the qp vacuum, and has a structure which can be described 
as a full quasiparticle shell. Notice that at this point the exact and full-qp 
excitation energies depart from each other. A state with four proton and 
four neutron quasiparticles has very large number-fluctuations. Spurious 
states become strongly mixed with physical states. In this way the resulting 
excitation energies average to zero, a limit which differs from the exact value, 
which is negative. 

The differences in the average qp number between full-qp and linearized 
approaches are larger in Figs. 6. a) and 6.b) than in Figs. 7. a) and 7.b). This 
result is a consequence of the dependence of some of the effective couplings 
of Eq. [ll], i.e: A3 ~ v? v — vi and A4 ~ — v%, on the number of particles. 
As the number of particles approaches the saturation value f2 = the 
effective couplings A 3 and A 4 will vanish. Thus, the full-qp and linearized 
solutions yield similar results. 



4 QRPA and RQRPA 

The QRPA hamiltonian Hqrpa can be obtained from the linearized version 
of Eq. (|TT|), by keeping only the bilinear-terms in the pair-creation and pair- 
annihilation operators. The pair-creation and pair-annihilation operators, 
and A, are, of course defined by coupled pairs of fermions. The commutation 
relations between these pseudo-boson operators include number-like quasi- 
particle operators in addition to unity. By taking the limit (2j + 1) — > 00 
|l"3|1 these extra terms vanish and the commutation relations between pairs 



of fermions can be treated like exact commutation relations between bosons. 
This is the well-known quasi-boson approximation and the QRPA hamilto- 
nian is the leading order hamiltonian which satisfies the quasi-boson approx- 
imation. If the pair-operators are replaced by quasi-bosons, the resulting 
Hamiltonian is given by 

Hqrpa = (2e + Ai) tfb + A 2 { 6+6+ + bb}. (22) 
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As said above, the quasi-bosons and b fulfill exactly the commutation 
relation [b, b'] = 1. 

At this point we can refer to pair of fermions ( A') or to quasi-bosons 
(6"!") without lost of generality, since we have not introduced a particular 
representation for the pair of fermions to boson mapping. 

The QRPA states are generated by the action of the one-phonon operator 
Oqrpa = XA^ — Y A on the correlated QRPA vacuum |0). The quasi-boson 
approximation assumes that (0|LA, -A']|0) = 1 and it leads to the normaliza- 
tion condition X 2 — Y 2 = 1. The QRPA matrix is just a 2 x 2 one, with 
sub-matrices Aqkva = 2e + Ai and Bqkva — 2A 2 . The corresponding eigen- 
value is given by Eqrpa = [(2e + A1) 2 — AX 2 ,} 1 / 2 . It becomes purely imaginary 
if2A 2 >2e + Ai. 

For this limit the backward-going amplitudes of the QRPA phonon-operator 
become dominant, thus invalidating the underlying assumption about the 
smallness of the quasi-boson vacuum-amplitudes. The QRPA excitation en- 
ergies, obtained with the above introduced Hamiltonian are shown in Figs. 
4 and 5. It can be seen that in the four cases displayed in these figures the 
collapse of the QRPA values occurs near the point where the exact excitation 
energies become negative. This is a very important result because it means 
that that the QRPA description of the dynamics given by the hamiltonian 
([I]) is able to reproduce exact results. At this point one can naturally ask the 
obvious question about the nature of the mechanism which produces such a 
collapse. The fact that the QRPA approximation is sensitive to it, together 
with the fact that the same behaviour is shown by the exact solution, rein- 
forces the idea about the onset of correlations which terminate the regime of 
validity of the pair-dominant picture. In order to identify such correlations 
we have calculated the expectation value of the number of quasi-fermions 
and bosons on the QRPA ground state. 

The average number of proton quasiparticles in the QRPA ground state, 
which in this case coincides with the average boson number, is given byQ 

(0|i\gO> = Y 2 . (23) 

Figs. 6. a), 6.b), 7.a) and 7.b) show the results corresponding to these 
occupation numbers. The QRPA results extend up to the value 4k/G ~ 1, 
where the QRPA collapses. The sudden increase of the average quasiparticle 



1 Notice that there is a factor 2 missprinted in Eq. (19) of ref. 13 
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number near the collapse of the QRPA is a clear evidence about a change in 
the structure of the QRPA ground state. 

In the renormalized QRPA the structure of the ground state is included 
explicitly ||, in the form 



|0) = Ne s \BCS) , S 



2(0|[AAt]|0) 



(24) 



where the quasi-boson approximation, at commutator's level, is not enforced 
explicitly. The renormalization procedure consists of retaining approximately 
the number of quasiparticle-like-terms of the commutators keeping them as 
a parameter to be determined, namely, by defining the RQRPA one-phonon 

'XA*-yA] /(0 |[A^ f ]|0) 1/2 1 0) (25) 



OrqrpaI®) 



and enforcing the condition Orqrpa\Q) = 0, which leads to the estimate c 
y j X for the parameter entering in the definition of the correlated vacuum 
After some algebra it is possible to show that (0|[A, A^O) = D 
[0, [TJ, and that 

2y 2 



-i _ 2y 2 D 

2j+l 



D =[ 1 + lTTj ' (26) 

The RQRPA submatrices are Ahquva — 2e + \\D and Bhqhp^ = 2\ 2 D. 
Since < D < 1, the presence of D multiplying both Ai and A 2 produces the 
reduction of the residual interaction which is needed to avoid the collapse of 
the QRPA equations 



Due to this fact the RQRPA energy Erqrpa is 
always real. Its value can be obtained by solving simultaneously the non- 
linear equations for Erqrpa, X, y and D, which in the general case will 



include all possible values of the multipolarity J [ITT 



RQRPA excitation energies are shown in Figs. 4 and 5. These results 
strongly resemble those of Fig. 1 of ref. JTT] and Fig. 2 of ref. M. The 
results corresponding to the QRPA (of ref.||) and to the RQRPA (of ref. 
[ |Tl|j ) are quite similar to those shown in Figures 4 and 5. However, the main 
finding of the present calculations is that the exact excitation energies are 
closer to the QRPA energies, rather than to the renormalized ones, instead. 
In exact calculations including the spin degrees of freedom a phase transition 
was found at the point where the QRPA collapses [p8 
present results. 



thus reinforcing the 
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The average number of quasiparticles in the RQRPA vacuum is given by 

(o|Ayo> = y 2 , (27) 

and it is shown in Fig. 6. a), 6.b), 7. a), and 7.b). It is fairly obvious, from 
these results, that the RQRPA ground state correlations double in all the 
cases those of the complete solutions of the linearized hamiltonian. This 
is clearly an overestimation, and it is probably one of the most notorious 
difficulties confronting the use of the RQRPA. 

It allows too much ground state correlations, and with them the particle 
number fluctuations are introducing spurious states which can dominate the 
low energy structure for large values of k. 

Near "collapse" the average number of quasiparticles given by the QRPA 
and the RQRPA are comparable. For the case of the QRPA the increase of 
the ground-state-correlations is determined by the change in the sign of the 
backward-going matrix relative to the forward-going one near collapse. From 
there on the QRPA cannot produce any physically acceptable result since one 
of the underlying conditions of the approximation, i.e: the positive definite 
character of either linear combination of the forward- and backward-going 
blocks of the QRPA matrix will not be fulfilled. This collapse is prevented 
in the RQRPA, by the use of the renormalization of the matrix elements, 
but the drawback of the approximation is the contribution coming from spu- 
rious states, which ought to be removed. Moreover, there are several other 
reasons to cast doubts on the consistency of the RQRPA. Among them, the 
mixing-up of orders in the wave functions, of the RQRPA phonons, is not 
accompanied by the the enlargement of the hamiltonian, to accommodate 
other correlations, like: i.e: the exchange terms of the QRPA matrix. If 
one performs such a calculation, by including exchange terms, the resulting 
values of the QRPA matrix terms are also " renormalized" , but this effect 
will depend upon the configurations. Also, the point of collapse is shifted 
to higher values of the coupling constant k but the effect is tipically of the 
order 1/fi, as compared to leading order terms. If terms others than unity 
are introduced in the commutators, then the hamiltonian has to be enlarged 
to account for the AB sort of terms of the initial hamiltonian, see eq.(ll), 
because they will contribute at the same order as the added number-type 
of terms introduced by the RQRPA procedure. Thus, the RQRPA proce- 
dure should be accompanied by a renormalization of the transition operators 
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and/or by the inclusion of scattering terms also in these operators. At this 
level, by going beyond the leading order QRPA approximation, more terms 
have to be added to the diagrams which represent the transition amplitudes. 
It has been done for a pure seniority model in ||29|| . This approach, for cor- 
relations between pairs of like-quasiparticles, is already cumbersome and it 
introduces an unmanageable number of contributions, both to the QRPA ma- 
trix as well as to the transition operators. For unlike-pairs of quasiparticles 
the situation can be even worse, since the complete algebra, which supports 
the expansions, cannot be defined in a subspace where scattering terms are 
replaced by c-numbers. More details about these aspects will be presented 
in a forthcoming publication. 



5 Double beta decay 



In this section we shall briefly discuss some of the consequences of the previ- 
ously presented approaches on the calculation of nuclear double-beta decay 
observables. In the following we shall focus our attention on the two-neutrino 
mode of the nuclear double-beta decay, since the matrix elements governing 
this decay mode are more sensitive to nuclear structure effects than the ones 
of the neutrinoless mode. As said in the introduction we shall consider only 
double-Fermi transitions. The nuclear matrix elements of the two-neutrino 
double-beta-decay M 2u can be written as: 

where |0j), |0a and )\0f) represent the initial, intermediate and final nuclear 
states participant of the virtual transitions entering in the allowed second- 
order weak-processes. The energies of the initial and intermediate states are 
Ei and E\, respectively. The energy released by the decay is represented by 
the quantity A. For the present calculations we have selected the value of 
A = 0.5 MeV, which is of the order of magnitude of the empirical values used 
in realistic calculations. The results for the matrix elements Mi v , obtained 
with the exact wave functions are shown, as a function of the ratio 4k/ G, 
in Figs. 8. a) and 8.b). These results have been obtained with the following 
set of parameters: j = 9/2, (M p = 2, J\f n = 8) —> (N p = 4, N n = 6) and 
X = and 0.04MeV r , respectively. The values shown in Figs. 9. a) and 9.b) 
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correspond to j = 19/2, (M p = 6, Af n = 14) -> (7V P = 8, M n = 12) and x = 
and 0.025MeV. 
Fig. 8,9 

For all cases the exact value of M 2v vanishes at the point 4k /G = 1. As 
mentioned above, this cancellation appears in the model due to the fact that 
for this value of k the isospin-symmetry is recovered and the ground states 
of the initial and final nuclei belong to different isospin multiplets, as it can 
be seen also from the results shown in Figs. 1 and 2. 

A similar mechanism, in the context of a solvable model possessing a 
SO (8) algebra including spin and isospin degrees of freedom was used a 
decade ago to show that the cancellation of the M 2v matrix elements for 
certain values of the particle-particle residual interaction was not an artifact 
of the QRPA description @. 

The results corresponding to the matrix elements M 2u , calculated with 
the different approximations discussed in the text are shown in Figs. 8 and 9, 
as a function of the coupling constant k. The values of M 2v are very similar 
to those found in realistic calculations [JXJ, [| f|, |11|, including its strong sup- 
pression for values of the coupling constant k near the value which produces 
the collapse of the QRPA description. Distinctively, the RQRPA results ex- 
tends to values of k passing the "critical" value. However, the validity of 
this result can be questioned because, as we have shown above, the RQRPA 
missed the vanishing of the excitation energy. The M 2v matrix elements, 
evaluated with the complete qp-hamiltonian (|il~l), is quite similar to that of 
the RQRPA up to point where it vanishes. From this point-on the results of 
both the full-qp and the RQRPA approximations are different. Both matrix 
elements change their sign at a value of k which is larger than the one cor- 
responding to the change of the sign of the matrix elements calculated with 
the exact wave function. The fact that the RQRPA results and the ones of 
the qp-approximation are similar, although these models differ drastically in 
the correlations which they actually include, suggest that a kind of balance is 
established between terms which are responsible for ground state correlations 
and those which produce the breaking of coherence in the wave functions. 
Obviously this mechanism must be related to the presence of scattering terms 
in the commutators as well as in the Hamiltonian. 
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6 Conclusions 



An exactly solvable model for the description of single- and double-beta- 
decay-processes of the Fermi-type was introduced. The model is equivalent 
to a complete shell model treatment in a single-j shell for the adopted hamil- 
tonian. It reproduces the main features of the results obtained in realistic 
calculations, with many shell and effective residual interaction, like those 
used in the literature to describe the microscopic structure of the nuclei in- 
volved in double beta decay processes. 

We have constructed the exact spectrum of the Hamiltonian and dis- 
cussed its properties. The results concerning the energy of the states be- 
longing to the exact solution of the model show that, in spite of its very 
schematic structure, the hamiltonian is able to qualitative reproduce the nu- 
clear mass parabola. The sequence of levels of the exact solution shows that 
the ground-state and the first-excited state, of the spectrum of double-even 
nuclei, approach a band-crossing situation for a critical value of the strength 
associated to attractive particle-particle interactions. At the crossing these 
states interchange their quantum numbers. This behaviour is connected with 
the description of " shape" transitions in similar theories, where the order pa- 
rameter is clearly associated with multipole deformations in r-space. In the 
present model the "deformation" mechanism is related with the breaking 
of the isospin symmetry and the space-rotation correspond to a rotation in 
isospin-space which preserves the third-component of the isospin. 

We have compared the exact values of the excitation energy and of the 
double-beta-decay matrix elements, for double-Fermi transitions, with those 
obtained by using the solutions of the approximate qp-hamiltonian, its lin- 
earized version and both the QRPA and RQRPA ones. 

It was shown that the collapse of the QRPA correlates with the presence 
of an exact-eigenvalue at zero energy. The structure of the RQRPA solutions 
has been discussed and it was found that though finite they are not free from 
spurious contributions. The role of scattering-terms was discussed and they 
were shown to be relevant in getting excitation energies closer to the exact 
values. However they are not enough to generate the correlations which are 
needed to produce the band-crossing, or negative excitation energies, as it 
was found in the exact solution for large values of the coupling constant k. 

In order to correlate the break-up of the QRPA approximation with the 
onset of strong fluctuations in the particle number we have calculated the 
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average number of quasiparticles in the different approximations discussed 
in the text. 

It was shown that the solutions of the complete qp-hamiltonian display a 
strong change in the structure of the ground state when the particle-particle 
strength increases. The qp-content of the ground state varies from a nearly 
zero-value to an almost full qp-occupancy. The particle number fluctuations 
associated with states with a large number of quasiparticles were mentioned 
as a possible source of spurious states. 

Double beta decay amplitudes were evaluated in the different formalisms. 
Their similitudes and differences were pointed out. 

As a conclusion the need of additional work, to clarify the meaning of the 
different approximations possed by the RQRPA, was pointed out. 
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Appendix A: the SO (5) algebra 



Following [g(J we introduce the operators 

j=o 



At(0) = 
fit 



Af=0 T 
J=0 



At, 



J=0 
M=0 



M=0 



«t <g> «t 



J=0 
M=0 



(29) 



which together with their hermitian conjugates and with the number and 
isospin operators 



N = N P + N n 



71 



N p -N n 



Ni 



^2at m .a im . i=p,n (30) 



are the ten generators of the SO (5) group. 

The hermitian conjugates of the pair-creation operators transform, under 
isospin-reversal like 

A(M) = (-1) M A(-M). (31) 

Their commutation relations are more easily expressed defining the new 
operators 



ffi = ?-n, 

E n = v/04t(l), 

= -VnAt(-i; 

£ 10 = v^At(O), 



#2 = 7z, _ 

= VnA{i), 

E- n = -VUA(-l), 

e_ w = VnA(o), 

l r 



(32) 



#0-1 = ~k T ■ 



The operators E a p are raising and lowering operators. When operating 
on an eigenstate of the weight operators H\ and H 2 they increase or decrease 
the eigenvalues of one or both by one unit. Their commutation relations are 



[Hi, H 2 ] = 0, [Hi, E a p] = aE a/3 , [H 2 , E a/3 ] = (3E a p, 
[E a p, E^a^p] = aHi + (3H 2 , , , 

f ± E a+a//3+l3 , if a + a' and P + p' = 0,1, -1 1 } 
I otherwise 



[E a p, E a >p>] 
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More explicitly 



[E n ,E_i ] = E i, [Exi,E -i] = —E 10 , [Eio,E-ix] = —E 01 , 
[Ex ,E-i-x] = E _x, [Exo,E 01 ] = —E X x, [E 10 , E _x] = E X -x, (34) 
[Ex-x, E-xo] = —Eq-x, [Ex-x,Eqx] = Exo 

and by hermitian conjugation of the above commutators one obtains 

ElL = E-a-p. (35) 



Appendix B: SO(5) representations 

The highest weights of the operators Hx , H2 define the irreducible rep- 
resentations (irrep) of the SO (5) algebra. For the present case we want the 
irrep which contains the state with zero quasiparticles as well as the state 
completely filled with quasi-proton and quasi-neutrons. The maximum num- 
ber (N max ) of quasiparticles allowed by the Pauli principle is 2Q, thus adding 
quasi-protons and quasi-neutrons one obtains N max = 40. This is state with 
the highest weight and it belongs to the irrep defined by (Hx = Q,H 2 = 0) 
ot N = 4Q, T — T z — 0. Acting with the generators ( |29| ) on this state it is 
possible to generate the set of all the states with even number of quasipar- 
ticles. This subspace suffices for all the calculations described in this work. 
For this reason we have adopted the irrep (Hx — £1, Hi = 0). 

In general it is necessary to specify four quantum numbers to completely 
define a state in a given irrep. But for the present case it turns out that the 
states can be defined by the quantum numbers N, T, T z . 

In the following we will construct explicitly the states \NTT Z = T); others 
states with T z ^ T are obtained by acting with the isospin lowering operator 
T~ on them. The states of this basis are defined by 



\N T = T z ) = 
0+ = E_ n , 
a = T z = T, 

N(a,b) = 2 b 



N(a,b)(Ooo) b (0 + r\N = 
Oqo = 1E_x\E-x-x + E. 
b = Q 



Ail T = 
iqE-xo, 



T z = 0), where 



T _ N 
2 4 ' 
(2Q+l-2by.(n-a-by.(2a+iy.(a+by. 

(2n+iy.(n-by.{a\) 2 b\(2a+2b+iy 



(36) 



1/2 
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Appendix C: SO (5) reduced matrix elements 



To diagonalize hamiltonian (11) in the (N,T,T Z ) basis, or hamiltonian 



([TJ) in the TV, T, T z basis, requires the use of the Wigner-Eckart theorem 

(N'T'T' z \O ntu \NTT z ) = (TT z ,tt z \T'T' z )(N'T'\\O nt \\NT) (37) 

where in the right hand side the symbol (.., ..|..) represents a Clebsch-Gordan 
coefficient and (.. 1 1 [..] 1 is a reduced matrix element. Explicit expressions for 
the reduced matrix elements are given below. The difference in the number 
of creation and annihilation operators in the tensor O is represented by n 
and in order to obtain non-zero matrix elements it must be equal to N' — N . 

We have used of the Wigner-Eckart theorem, the commutation relations 
given in the Appendix A and the explicit form of the states with T = T z 
shown in the Appendix B to calculate the reduced matrix elements of the 
different operators which are relevant in our problem. Some of these SO (5)- 
reduced matrix elements are listed here. Additional matrix elements can be 
deduced from them by using 



2T' 4- 1 

(NT\\{O nt y\\N'T') = \ (N' T'\\O nt \\NT) . (38) 

y LI -\- 1 

The relevant reduced matrix elements are 



(iVT + 2||[Atir= 2 ||iVr) = ^ 

(2tt - T - N/2)(T + N/2 + 3)(-T + iV/2)(2ft + T - N/2 + 3)(T + 1){T + 2) 



1/2 



(2T + 3)(2T + 5) 



(NTW^Af^WNT) = — [(2fl - N/2 + 3)N/2 - T(T + 1)] , 



(NT\\[A^AY= 2 \\NT) = ^—[{NT = T z \A\l)A{l)\NT = T Z )+ 

(NT = T z \A j (-l)A(-l)\N T = T Z ) — 2(N T = ^1^(0)^(0)1^ T = T z ) 
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where 



(N T = T Z \A*(1)A(1)\N T = T z ) 



1 

ft 



-ft + T + N/2 + 



(2ft - T - JV/2)(T + iV/2 + 3)(T + 1) 
2(2T + 3) 



(N T = T Z \A\-1)A(-1)\N T = T z ) 



1 

ft 



(2ft — T — JV/2 + 3)(-T + N/2)(T + 1) 
2(2T + 3) 



(N T = T Z \A* (0)A(0)\N T = T z ) = - [-ft + JV/2- 



(2ft — T — iV/2)(T + iV/2 + 3)ft 
(2ft + T — iV/2 + 1)(-T + iV/2 + 2 



■(iV + 4T = ^1^(0)^(0)1^ + 4 T 



The largest value that N can take is 4ft — 2T. In this case 

(N = 4ft - 2T T = ^1^(0)^(0)1^ = 4ft - 2T T = T z ) = 1 - T/ft 

The above reduced matrix elements are enough to deal with hamilto- 
nian (]T|), which conserves particle number. Working with hamiltonian (|TT|) 
requires many other reduced matrix elements, like the following matrix ele- 
ments 



(N + AT\\[A^A^ t=0 \\NT) 



-1 



[(T + N/2 + 3)(-T + N/2 + 2) (2ft - T - iV/2)(2ft + T-N/2 + 1)] 



1/2 



(iV + 4 T - 2\\[A ] A^} t=2 \\NT) 



2ft 



(-T+N/2+4)(-T+N/2+2)(T-l)T(2n+T-N/2-l)(2Q+T-N/2+l) ] 1 / 2 
(2T-l)(2T-3) 
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(N + AT + 2\\[A ] A ] ] t=2 \\NT) = — 

' (T+N/2+3)(T+N/2+5)(T+l)(T+2)(2Q-T-N/2)(2n-T-N/2-2) 1 1 / 2 



(2T+3)(2T+5) 



(iV + 2 r-lH^HJVT) 



T(-T + iV/2 + 2)(2fi + T - AT/2 + 1) 
2fi(2T - 1) 



(T + 1)(T + iV/2 + 3) (2fi — T — JV/2) 
2fi(2T + 3) 



nl/2 



(AT + 2 T + l\\A*\\NT) 



These matrix elements, together with those associated with the isospin 
raising and lowering operators 

T+ = , T~ = -V2£LB ] , 

(NTT Z + l\B\N T T z ) = — [(T + T z + 1)(T - T 2 )] 1 / 2 /v / 2Q , 
(N T T z — 1\B*\N T T z ) = — [(T - T* + 1)(T + T^p^/v^ , 

are all the elements which are needed to diagonalize the hamiltonian (|TT|) 
and to calculate the matrix elements of the transition operators. 
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Figure Captions 



Figure l.a (l,b): + states of different isotopes are shown for j = 9/2, 
4k/ G = 1 and X — 0.(0.05) MeV, in an energy vs. Z plot. States are labeled 
by (T,T Z ). The lowest energy state of each nucleus is shown by a thick-line. 

Figure 2. a (2.b): The same as Fig 1 for j = 19/2, 4k/G = 1 and X = 
0. (0.025)Me\/. 

Figure 3. a): Energy of the ground state 0+ (full line) and first excited 
state 0i (dotted line), as a function of the ratio 4k/ G, for j = 9/2, M n = 
6, Af p — 4, x — 0- Figure. 3.b) shows the same quantities for j = 19/2, Af n = 
12, J\f p = 8. 

Figure 4. a (4.b): Excitation energy E exc of the lowest + state in the 
odd-odd intermediate nucleus (N n = 7,J\f p = 3) with respect to the parent 
even-even nucleus (Af n = 8,J\f p = 2) against 4k/ G for j = 9/2, x = (0.04). 
Exact results are shown as thin-full-lines while those of the qp-hamiltonian 
are shown as small-dotted-lines. Results corresponding to the linearized qp- 
hamiltonian are shown as full-thick-lines and the results obtained with the 
QRPA and RQRPA methods as large-dotted- and dashed-lines, respectively. 

Figure 5. a (5b): The same as Figure 4. a) but for the excitation energy E exc 
of the lowest + state in the odd-odd intermediate nucleus (Af n = 13, Ap = 7) 
with respect to the parent even-even nucleus (Af n = 14, J\f p = 6) , for j = 
19/2, x = (0.025) 

Figure 6. a (6.b): Average number of proton-quasiparticles in the ground 
state of the even-even nucleus with M p = 6, J\f n = 14 as function of 4k/ G, for 
j = 19/2, x = (0.025)MeV . Results corresponding to the qp-hamiltonian 
are shown as dashed-lines. The ones corresponding to the linearized qp- 
hamiltonian are shown as large-dotted lines and those of the QRPA and 
RQRPA methods as full-lines and small-dotted-lines, respectively. 

Figure 7. a (7b): The same as Figure 6 for the number of proton-quasiparticles 
for M p = 8, N n = 12, j = 19/2, x = (0.025)Me1/. 

Figure. 8a (8b): Matrix elements M 2u , for the double- Fermi two-neutrino 
double-beta decay mode, as functions of the ratio 4k/G for j = 9/2, (Af p = 
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2, M n = 8) -> (jV p = 4, JV n = 6) and x = (0.04)MeV. Exact results are 
indicated by thin-full-lines. The results obtained with the qp-hamiltonian 
are shown as small-dotted-lines and the results of the QRPA and RQRPA 
methods as dashed-lines and large-dotted-lines, respectively. 

Figure 9. a (9.b): The same as Figure 8), i.e: the matrix elements M 2u , for 
j = 19/2, (Up = 6, N n = 14) -> (N p = 8, N n = 12) and x = (0.025)MeV 
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Figure 7 
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Figure 8 
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Figure 9 
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